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ABSTRACT 


\ 

^Tbe  Wilson  Wind-Generated  Source  Level  Densities 
and  the  ASTRAL  propagation  loss  model  are  combined  to  form 
a  wind-generated  noise  model.  Tbe  model  treats  storms  as 
a  geographically  distributed  source  and  is  capable  of 
modeling  both  distant  and  local  storm  noise.  Model  es¬ 
timates  are  compared  with  omnidirectional  measurements 
from  the  northeastern  Pacific  Ocean  at  5  frequencies  and 
3  receiver  depths.  The  results  are  very  good  and  indicate 
the  necessity  of  modeling  both  shipping  and  wind-generated 
noise  in  comparisons  with  measured  data. 

/ 


IV 


1.0 


INTRODUCTION 


An  ambient  noise  model  has  three  components. 
First,  a  geographical  distribution  of  the  source  o f  noise 
(e.g.,  shipping,  wind,  etc.)  oust  be  specified  along  with 
the  geographical  location  and  depth  of  the  receiver. 
Second,  the  source  level,  usually  specified  at  some  ref¬ 
erence  distance  from  the  source,  and  source  depth  and 
directivity  pattern  must  be  specified.  Third,  an  accurate 
propagation  loss  model  must  be  applied  between  the  source 
distribution  and  receiver  in  order  to  estimate  ambient  noise 
levels  at  the  receiver  site.  The  now  familiar  Wenz*  curves 
are  not  an  ambient  noise  model  in  this  sense.  They  are 
noise  levels  representative  of  deep  water  ocean  basins. 
Recent  acoustic  and  oceanographic  at-sea  measurements 
made  under  the  sponsorship  of  the  Naval  Oceanographic 
Research  and  Development  Activity  have  shown  the  highly 
variable  nature  of  ambient  noise  level  as  a  function  of  time 
and  source/receiver  location.  One  set  of  ambient  noise 
level  curves,  while  a  valuable  reference,  does  not  provide 
the  real  time  ambient  noise  input  that  is  needed  for  the 
optimization  of  performance  of  current  tactical  and  surveil¬ 
lance  sonar  systems.  The  wind-noise  model  described  in  this 
paper  does  provide  a  dynamic,  real  time  estimation  of  noise 
due  to  storms  at  sea. 

In  the  author’s  opinion,  there  are  three  signif¬ 
icant  applications  of  this  type  of  wind-ooise  model.  First, 
at  any  given  time  the  model  cao  identify  high  noise  areas 
and  directions  where  search  for  signals  of  interest  will  be 
unsuccesful  with  high  probability.  Thus,  sonar  system 
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resources  can  be  allocated  accordingly.  Conversely,  the 
model  can  estimate  where  and  in  what  directions  the  noise 
level  due  to  storms  is  low  so  that  mobile  sonar  systems  can 
be  placed  in  the  optimum  location,  depth,  and  heading  to 
maximize  probability  of  detection  of  signals  of  interest. 
Second,  the  model  can  forecast  noise  levels  up  to  3  days  in 
the  future  if  the  storm  path  and  intensity  can  be  forecasted 
accurately.  This  can  be  used  to  estimate  what  the  optimum 
sonar  system  location  and  heading  will  be  in  the  future,  or 
what  sonar  system  resources  will  be  necessary  to  search, 
localize,  and  maintain  contact  with  a  signal  of  interest. 
Third,  if  storm  noise  can  be  accurately  modeled,  there  is  a 
possibility  of  subtracting  the  noise  from  sonar  system 
measurements  to  improve  system  performance.  This  as  yet, 
has  not  been  demonstrated  and  will  be  the  subject  of  further 
research. 


In  the  frequency  range  of  interest  for  this 
model,  10  to  1000  Hz,  shipping  and  wind/storm  noise  are  the 
major  components  of  the  ambient  noise  field.  As  described 
above,  wind  noise  will  be  modeled  from  real  time  wind  speed 
distributions  over  the  ocean  surface.  In  order  to  evaluate 
the  wind  noise  model  with  measured  data,  shipping  needs  to 
be  modeled  also  because  both  wind  and  shipping  contribute  to 
the  measured  data.  Historical  or  RMS2  shipping  densities 
will  be  used  in  the  model,  because  real  time  shipping 
densities  are  not  currently  available  as  model  inputs.  The 
model  output  will  consist  of  a  wind  noise  estimate  based  on 
dynamic  wind  speed  distributions  and  a  shipping  noise 
estimate  based  on  historical  shipping  densities  that  are 
constant  in  time. 
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The  wind-generated  noise  source  characteristics 
are  described  in  the  following  sections.  Details  of  the 
wind  noise  model  development  are  given  in  Section  3.0. 
Model/data  comparisons  are  described  in  Section  4.0 
for  data  collected  in  the  northeastern  Pacific  Ocean  in 
September  1973. 

2.0  WIND-GENERATED  SOURCE  LEVEL  DENSITY 

The  source  level  densities  of  wind-generated 
noise  are  shown  as  a  function  of  frequency  and  wind  speed  in 
Figure  1.  These  curves  are  the  result  of  analytical  studies 
of  the  physical  mechanisms  of  wind-generated  noise  from  5  to 
1000  Hz  performed  by  the  author. 3 Deep  ocean  measure¬ 
ments  by  Wittenborn  and  Shooter**  and  liorris®  compare 
favorably  with  the  results  shown  in  Figure  1. 

In  deriving  the  source  level  density  shown  in 
Figure  1,  great  care  was  taken  to  address  the  difference 
between  noise  level  and  source  level.  In  order  to  validate 
wind-generated  noise  source  levels  with  measured  noise 
levels,  only  a  very  select  set  of  measured  data  can  be 
considered.  In  the  Appendix,  it  is  shown^  that  source 
level  and  noise  level  are  related  by  a  constant  factor  only 
when  the  noise  measurements  are  taken  well  below  critical 
depth  near  a  region  with  a  poor  reflecting  bottom  where  only 
the  noise  from  the  near  vertical  sea  surface  above  the 
hydrophone  is  measured.  In  this  case,  the  relationship 
between  source  level  and  noise  level  for  a  dipole  radiation 
pattern  is  given  by 

NL(f,U)  *  *SL(f,D)  (1) 
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where  f  is  frequency  in  Hz  and  U  is  wind  speed  in  knots. 

In  most  measured  noise  data,  where  the  hydrophone 
is  above  the  critical  depth,  noise  from  distant  sources 
becomes  important  and  the  frequency  dependence  of  the 
transmission  loss  and  bottom  interaction  make  it  impos¬ 
sible  to  relate  noise  level  and  source  level  by  the  simple 
relation  given  in  Equation  (1). 

Acoustic  Data  Capsule  (ACODAC)  data  from  the 
northeastern  Pacific  analyzed  by  Wittenborn  and  Shooter5 
and  FLIP  data  from  the  same  area  analyzed  by  Morris6  were 
used  to  validate  the  analytical  results  shown  in  Figure  1. 
The  measured  data  and  analytically  derived  levels  compare 
favorably.  However,  more  of  these  near  bottom  measure¬ 
ments  are  needed  for  a  wider  range  of  wind  speeds  and 
frequencies,  since  not  all  frequencies  and  wind  speed 
intervals  were  covered  by  the  measured  data. 

Figure  1  shows  noise  level  values  as  a  function 
of  frequency  and  wind  speed  and  the  corresponding  source 
level  values  can  be  derived  by  using  Equation  (1)  (i.e., 
divide  the  noise  level  value  shown  in  Figure  1  by  it  ) . 
Table  I  lists  the  source  level  values  derived  from  Figure  1 
at  frequencies  selected  for  input  to  the  model.  Linear 
interpolation  between  the  values  shown  in  Table  I  is  used  in 
the  model  for  frequencies  and  wind  speeds  not  listed  in 
Table  I. 


Table  I 

WIND-GENERATED  SOURCE  LEVEL  DENSITIES  (dB/y Pa/Hz1/2) 
FOR  DIPOLE  RADIATION  PATTERN 
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2.1 


CHARACTERISTICS  OF  THE  WIND-GENERATED  NOISE 
SOURCE 


Analyses  of  the  physical  mechanisms  of  wind¬ 
generated  noise  indicate  that  turbulene  pressure  fluctua¬ 
tions  in  the  atmosphere  near  the  ocean  surface  was  the 
dominant  physical  mechanism  from  10  to  80  Hz,  while  wave- 
wave  interaction  was  the  dominant  physical  mechanism  from  5 
to  10  Hz  and  the  impact  of  ocean  spray  streaks  and  whitecaps 
with  the  ocean  surface  was  the  dominant  physical  mechanism 
at  higher  frequencies  from  80  to  1000  Hz.  All  of  the 
above  pnysical  mechanisms  had  a  dipole  radiation  pattern 
associated  with  them.  The  exact  frequency  interval  over 
which  a  particular  mechanism  is  dominant  is  highly  wind 
speed  dependent  and  the  above  intervals  are  only  rough 
approximations. 

The  wind  speed  dependence  of  the  wind-generated 
source  levels  was  significantly  different  at  low  and  high 
frequencies.  From  5  to  80  Hz,  the  wind  speed  dependence  of 
the  analytically  derived  source  levels  and  relevant  measured 
noise  levels®  varied  as  wind  speed  to  the  3.65  power, 
U®*65.  This  result  was  also  confirmed  by  Perrone's® 
measurement  in  the  Atlantic  Ocean.  At  frequencies  above 
200  Hz,  the  wind  speed  dependencies  of  the  analytically 
derived  source  levels  and  relevant  measured  noise  levels'* 
were  directly  proportional  to  the  so-called  "whitecap 
index,”  R(U).  The  whitecap  index  is  defined  as  the  percent 
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of  the  ocean  surface  covered  by  whitecaps  and  streaks  and 
varies  as  a  function  of  wind  speed  (in  knots)  as 


R(U)  = 


'•  .0 

U3  -  U2  +JL_-1.5 
1749.60  81.00  4.32 

-  R(30)(U/30)L5 


0  kts  <U<9  kts 
9  kts  <U<30  kts 

.  <30  kts  (  2  ) 


Since  the  onset  of  whitecaps  occurs  for  wind- 
speeds  of  9  knots,  R(U)  was  taken  to  be  zero  below  9  knots. 
For  the  wind  speed  interval  between  9  and  30  knots,  R(0) 
was  empirically  fit  to  experimental  measurements,9  and 
above  30  '-nots  the  saturation  of  R(U)  was  taken  into  account 
using  Perrone's  data.®  This  saturation  effect  was  not 
considered  i:  the  author's  latest  published  analyses.4 


The  frequency  dependence  of  the  wind- generated 
source  level  was  relatively  flat  from  10  to  1000  Hz  where 
turbulent  pressure  fluctuations  in  the  atmosphere  and  the 
impact  of  ocean  spray  and  streaks  were  dominant  physical 
mechanisms.  A  negative  slope  is  predicted  from  5  to  10  Bz 
where  wave-wave  interaction  is  the  dominant  physical  mech¬ 
anism.  The  measured  data^>6  show  a  relatively  flat 
spectrum  above  10  Hz,  but  there  is  no  data  available  to 
validate  the  negative  slope  predicted  from  5  to  10  Hz. 


The  source  level  densities  presented  in  Fivtre  1 
are  only  a  first  estimate  of  the  wind-generated  source 
levels  that  are  consistent  with  the  appropriate  measured 
data.  Many  further  analytical  developments  and  appropriate 
measurements  are  necessary  to  refine  and  completely  validate 


wind-generated  source  levels.  The  source  levels  shown  in 
Figure  1  are  substantially  different  than  the  Wenz  curves1 
because  the  former  is  an  estimation  of  source  level  while 
the  latter  is  an  estimation  of  the  noise  level.  The  Wilson 
source  level  curves  are  intended  for  use  with  an  appropriate 
transmission  loss  model -to  estimate  ambient  noise  levels. 

3.0  DEVELOPMENT  OF  WIND-GENERATED  NOISE  MODEL 

The  wind-generated  source  level  densities 
presented  in  the  previous  section  need  to  be  combined  with 
an  appropriate  transmission  loss  model  in  order  to  estimate 
ambient  noise  levels.  ASTRAL10  was  selected  because  it  is 
currently  being  used  in  the  model  to  estimate  shipping 
noise,  and  because  ASTRAL's  range  averaging  property  is  very 
compatible  with  a  source  such  as  wind,  which  is  distributed 
over  large  ocean  areas. 

In  any  ambient  noise  model  that  includes  both 
wind  and  shipping,  it  is  efficient  to  treat  both  noise 
sources  alike  when  applying  the  transmission  loss  model. 
To  accomplish  this  objective,  the  concept  of  a  "Wind  Equiv¬ 
alent  Ship"  (WES)  was  developed.  Since  the  model  computes 
ship  noise  from  shipping  densities  in  1°  x  1°  Latitude/ 
Longitude  cells,  wind  source  level  densities  were  calculated 
for  the  same  1°  x  1°  cells.  This  was  done  simply  by 
multiplying  the  area  (Aiyd2)  of  the  itb  1°  x  1°  cell 

A*  »  8.25  x  1011  (cos  L£  -  cos  Lu)  (3) 
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(where  Lj,  and  Lu  are  the  latitudes  in  the  1°  x  1°  cell 
nearest  to  and  farthest  from  the  equator,  respectively)  by 
the  source  level  density  from  Table  I,  SL(Ui,f),  for  the 
frequency  of  interest  and  for  the  average  wind  speed,  Ui, 
over  the  i*k  cell: 

SL(Ui,f)  -  SL(Ui,f)  Ai  .  (4) 

The  source  level,  SL(Ui,f),  is  equivalent  to  a  ship 
of  the  same  source  level  located  at  the  center  of  the 
1°  x  1°  cell.  One  important  difference  between  a  ship 
and  the  WES  of  the  same  source  level  remains:  source  depth. 
Wind  is  a  near  surface  source  while  ships  radiate  noise  from 
a  point  well  below  the  surface  (a  20-ft  ship  source  depth  is 
assumed  in  the  model).  However,  since  the  wind  source 
directivity  is  believed  to  be  dipole  in  nature,  its  source 
depth  can  be  lowered  to  depth,  z,  below  the  ocean  surface  by 


SL(Oiff,z)  * 


$L(D1f ,z-0) 
4k2z2 


(5) 


where  k  is  the  acoustic  wave  number.  As  shown  in  the 
appendix,  the  term  4k2z2  accounts  for  surface  image 
interference  for  kz  sins  <<  1. 


To  compute  the  ambient  noise  from  the  ith  cell 
as  measured  by  a  receiver  at  depth  d,  the  transmission  loss 
is  computed  between  the  source  coordinates  (source  location 


and  depth)  and  receiver  coordinates  (receiver,  location  and 
depth)  at  the  frequency  of  interest.  The  noise  from  the 
jth  Ceil  is 

Nt  -  TLi  !E(Ui,f,z)  (6) 

The  shipping  noise  is  computed  similarly  and  added  to  the 
wind  noise  and  only  one  transmisison  loss  computation  is 
made  for  the  itb  cell.  By  including  only  cells  from  a 
given  azimuthal  sector,  the  model  can  estimate  the  noise 
arriving  from  only  that  sector.  By  making  the  sector 
small,  the  noise  horizontal  directionality  can  be  estimated 
from  all  sectors  from  0°T  to  360°T.  If  the  measuring 
apparatus  is  an  array,  the  beam  pattern  can  be  convolved 
with  the  sector  noise  estimates  to  estimate  beam  noise.  Or, 
if  all  the  sectors  are  added  together,  the  omnidirectional 
noise  can  be  estimated.  Wind  noise  and  shipping  noise  can 
also  be  computed  separately.  Cells  that  are  split  between 
two  sections  have  their  wind  ard  shipping  source  levels 
split  between  the  sectors  proportional  to  the  area  split. 

The  above  procedure  is  currently  followed  for 
cells  greater  than  50  nm  from  the  measurement  site.  The 
resolution  of  the  transmission  loss  model  is  not  good  enough 
for  the  nearby  cells,  and  tbe  local  wind  noise  levels  are 
taken  from  Figure  1  directly  and  added  to  the  distant 
noise. 


4.0  MODEL/ DATA  COMPARISON 

ACODAC  data  collected  over  a  10  day  period  in 
September  1973  from  3  hydrophones  (696m,  4055m,  and  5521m) 
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located  near  39°N  Latltude/143°W  Longitude  are  used  in 
the  model/data  comparison.  Time  series  and  10  day  averages 
of  spectrum  levels  were  computed  at  12.5,  25,  50,  100,  160 
and  250  Hz.  Sound  channel  axis  depth,  critical  depth,  and 
water  depth  are  655m,  3860m,  and  5555m  respectively.  The 
average  local  wind  speed -for  the  10  day  period  was  20  knots 
and  there  were  no  major  dist^t  storms  during  this  period. 

Figure  2  shows  a  comparison  of  model  estimates 
with  10  day  average  spectrum  levels  at  five  frequencies  and 
three  receiver  depths  throughout  the  water  column.  Shipping 
noise  was  dominant  for  the  upper  two  receivers  (696m  and 
4055m)  below  100  Hz.  Wind  noise  was  dominant  for  these 
receivers  at  160  Hz  and  250  Hz.  The  model  estimates  for  the 
near  bottom  receiver  (5521m)  are  plotted  in  Figure  2  for 
wind  only,  shipping  only,  and  total  noise.  Wind  noise  is 
dominant  throughout  the  spectrum,  especially  above  100  Hz. 
Notice  the  flat  spectra  characteristic  of  wind  noise  in  both 
the  model  estimates  and  measured  data.  The  model/data 
comparisons  shown  in  Figure  2  are  very  good,  in  the  author’s 
opinion,  and  show  that  the  model  does  well  for  long  time 
averages.  How  the  model  estimates  compare  with  the  noise 
time  series  is  discussed  next. 

The  time  series  for  the  middle  hydrophone 
(4055m)  at  250  Hz  is  plotted  every  12  hours  in  Figure  3. 
If  a  nearby  ship  or  ship  tow  source  operations  were 
present  a  range  of  values  Is  plotted.  For  example,  on 
18  September  1973  a  nearby  ship  was  seen  in  the  250  Hz  time 
series  and  the  noise  level  increased  from  66  dB/  vPa/Hz*/2 
to  69  dB/PPa/Hz1/2  during  the  passage  of  this  ship.  This 
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Figure  2.  Model /Data  Comparison  -  10  Day  Averages 
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range  of  values  is  plotted  on  Figure  3.  The  local  wind 
spaed  for  each  12  hour  period  is  specified  along  the  time 
axis  of  Figure  3  and  the  model  wind  noise  estimates  are 
plotted  throughout  this  period.  The  model  estimate  of 
historical  shipping  noise  does  not  change  during  this  period 
since  the  model  shipping  -densities  and  model  transmission 
loss  is  constant.  The  model  estimates  of  total  noise  (wind 
plus  shipping)  is  plotted  and  compares  very  favorably  with 
the  measured  data  for  this  10  day  period. 

5.0  CONCLUSIONS 

This  preliminary  evaluation  of  the  wind  noise 
model  (and,  thus,  the  source  level  densities  given  in  Fig¬ 
ure  1  and  Table  I)  shows  very  good  agreement  be' .veen  model 
and  data,  in  the  author's  opinion.  However,  a  much  more 
extensive  evaluation  for  different  receiver  locations  and 
larger  storms  needs  to  be  performed.  This  preliminary 
evaluation  is  only  intended  to  be  a  start  in  the  right 
direction  in  modeling  storm  noise. 
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Appendix 

RELATIONSHIP  BETWEEN  NOISE  LEVEL  AND  SOURCE  LEVEL 

In  deriving  the  source  level  densities  in  Fig¬ 
ure  1,  great  care  was-  taken  to  address  the  difference 
between  noise  level  and  source  level.  .  Xn  order  to  validate 
wind-generated  noise  source  levels  with  measured  noise 
levels,  only  a  very  select  set  of  measured  data  can  be 
considered.  It  is  shown7  below  that  source  level  and 
noise  level  are  related  by  a  constant  factor  only  when  the 
noise  measurements  are  taken  well  below  critical  depth  near 
a  region  with  a  poor  reflecting  bottom  where  only  the  noise 
from  the  near  vertical  sea  surface  above  the  hydrophone  is 
measured.  In  this  case,  the  relationship  between  source 
level  and  noise  level  for  a  dipole  radiation  pattern  is 
given  by 


NL(f,U)  -  *SL(f,U)  (A-l) 

Assuming  that  only  direct  paths  are  important  for 
sea  surface  above  a  very  deep  receiver,  the  noise  level,  NL, 
and  source  level,  SL,  are  related  by 

NL  *  I  3L(8 0)  .  TL(80,6r)  dA  (A-2) 

Jocean 

Surface 


A-l 


where  the  transmission  loss  is  given  by 


C08  0q 

TL(0o,6r)  -  dR 

R  sln9r 


(A-3) 


and  where  the  ocean  surface  incremental  area  is  given  by 


dA  *  R  dR  d  4> 


(A-4) 


In  Equations  (A-2),  (A-3)  and  (A-4),  the  variables 
are  defined  as  follows: 

0q  *  surface  source  angle 
0r  =  receiver  arrival  angle 
R  ®  horizontal  source/receiver  range 
4>  *  azimuth  angle 

A  »a  ocean  surface  area 

Substituting  Equations  (A-3)  and  (A-4)  into  Equa¬ 
tion  (A-l>  and  integrating  over  azimuth  yields 


SL(eJ  cos(6J  d8n 
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with  Cr  and  Cq  being  the  sound  speed  at  the  receiver  and 
surface  respectively.  For  the  purpose  of  the  analysis,  a 
homogenous  (Cr  *  Cq)  medium  will  be  assumed  and  for  a 
•’sine  power"  source,  SL(6  o)  ®  SL(sinm  0o)»  Equation  (A-6) 

NL  *  *m~  '  (4-7) 

For  a  dipole  source  (m®2),  Equation  (A-7)  becomes 

SL  -  —  (A-8) 

TT 

Equation  (A-8)  is  the  connection  between  the  noise  levels 
plotted  in  Figure  1  and  the  source  levels  listed  in  Table  I. 

If  computing  efficiency  was  nu;  of  concern,  Equa¬ 
tion  (A-8)  C'uld  be  used  as  the  final  relationship  between 
source  level  and  noise  level.  However,  ASTRAL  transmission 
loss  is  currently  computed  for  shipping  in  the  model  for  a 
source  depth  of  20  feet.  In  order  to  avoid  another  computa¬ 
tion  cf  transmission  loss  for  the  wind  surface  (za»0)  source, 
the  following  source  depth  relation  is  used 

TL(z)  TL(z*0)  4  sin2(kz  sine  q)  (A-9) 


and 


for  kz  sine  <<  I 
(sin2  8  0)  TL(z=0)  » 

4k*z2 


( A-10) 


The  approximation  m  going  from  Equation  (A-9)  to  (A-10)  is 
valid  for  long  distance  (small  angles,  @  q) ,  low  frequency 
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calculations.  For  example,  for  f  =  150  Hz  (or  X  «  32.5  ft 
and  k  =  0.193  ft”*)  and  for  z  =  20  ft,  the  restriction  on 
source  angle  is 


8  0  <<  15°,  ( A- 1 1 ) 

which  is  a  good  approximation  for  long  range  propagation. 
Below  150  Hz,  the  approximation  gets  better,  and  Equa¬ 
tion  (A-10)  shows  that  a  surface  dipole  can  be  approximated 
by  a  monopole  source  at  depth  z  whenever  kz  sin  0  $  <<  1. 
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